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We present results of test-particle simulations on both the first- and the second-order Fermi acceleration for 
relativistic parallel shock waves. Our studies suggest that the role of the second-order mechanism in the 
turbulent downstream of a relativistic shock may have been underestimated in the past, and that the stochastic 
mechanism may have significant effects on the form of the particle spectra and its time evolution. Figures are 
reproduced from Virtanen & Vainio 2005 Hi by the permission of the A AS. 


1. INTRODUCTION 

Particle acceleration in relativistic shocks is typi¬ 
cally considered to be due to either the first- or the 
second-order Fermi acceleration (stochastic accelera¬ 
tion). The former of these acts at the shock front 
where charged particles gain energy by subsequently 
scattering on different sides of the shock front, whereas 
for the latter case the velocity difference between the 
scattering centers is provided by the turbulence alone. 
The first-order mechanism is well known to produce 
power-law particle energy spectrum N(E ) oc E~ a 
with spectral index <j depending on compression ratio 
r of the shock as cr = (r + 2)/(r — 1) for nonrelativistic 
speeds [e.g.,|lj, and approaching to value a ~ 2.2 as 
the speed increases to ultrarelativistic, almost regard¬ 
less of the other shock conditions [e.g.,|2, J, [j, Ej. 

However, problems arise with flatter spectra (cr < 
2); although the first-order mechanism can produce 
harder spectra with cr —> 1 depending, for instance, on 
the injection modellS] or the scattering center com¬ 
pression ratio [e.g. GlilHj it is n °t able to produce 
spectral indices flatter than cr = 1. 

Stochastic acceleration, on the other hand, has been 
known to be present in the turbulent downstream of 
shocks, but mostly because it works on much longer 
timescales than the almost instantaneous first-order 
mechanism [e.g., Him it has been been neglected 
for most of the cases (note, however, e.g. Q2)- How¬ 
ever, for non-thermal particle distributions radiating 
in astrophysical objects the bulk of radiation is emit¬ 
ted by the particles that have already left the shock 
front towards the downstream. Thus, the second or¬ 
der mechanism has, indeed, more time available to 
for acceleration than the first-order process. So while 
the neglecting of the stochastic process in calculation 
of the accelerated particle spectra right at the shock 
front could be justified, it is not possible to neglect its 
effect on the spectrum in shocks in general, especially 
when basically all of the astrophysical sources related 
to relativistic shocks, remain still spatially unresolved. 

Here we study the possibility of the stochastic pro¬ 
cess to have visible effects on the particle energy spec¬ 
trum in parallel relativistic shock waves in two differ¬ 


ent cases: (i) particles injected at the shock front and 
accelerated further by the first order mechanism and 
(ii) particles drawn from the heated (but not shock 
accelerated) particle population of the downstream re¬ 
gion of the shock. We apply numerical Monte Carlo 
test-particle simulations and focus on shocks that, in 
addition to being parallel, have small-to-intermediate 
Alfvenic Mach number. Low Mach numbers for rel¬ 
ativistic shocks could prevail in magnetically domi¬ 
nated jets that are lighter than their surroundings, 
e.g. in pair-plasma jets. 

2. MODEL 

The applied model, together with its numerical im¬ 
plementation, is described in detail in ref. 0 , and 
here only the essential parts of the model are ex¬ 
plained. 

We use the shock rest frame as our basic coordinate 
system. In the case of parallel shock the flow direc¬ 
tion - as well as the large scale magnetic field - is 
perpendicular to the shock normal in this frame. In 
this frame the shock lies at the origin, negative values 
of location x mean the upstream and positive values 
the downstream. The plasma frame, i.e. the frame 
where the bulk of the plasma is at rest, is moving 
with the local flow speed V with respect to the shock 
frame. Waves propagate at phase speed V# in the 
plasma frame in directions both parallel and antipar¬ 
allel to the flow, so in the shock frame waves move at 
speed (V + y>)/(l + HV^/c 2 ) (c is the speed of light), 
and the frame moving at this speed (i.e. flowing with 
the waves) is called the wave frame or the rest frame 
of the scattering centers. If the scattering centers are 
taken to be fluctuations frozen-in to the plasma then 
the speed of the waves with respect to the underlying 
flow is Vtf, = 0 and the plasma frame is also the rest 
frame of the scattering centers. Here we assume the 
waves to be Alfven waves, and the phase speed to be 
the Alfven speed 
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Figure 1: Ratio of the amplified wave intensities as 
function of Alfvenic Mach number for three proper 
Alfven speeds. Here q = 5/3. 0 


where Bq is the large scale magnetic field, and h and 
n the specific enthalpy and the number density, all in 
the local plasma frame. 

We set the proper speed of the shock to be u\ = 
TiVi = 10c (corresponding to shock speed V\ ~ 
0.995 c), and model the shock transition with a very 
thin hyperbolic tangent profile 0 with thickness 
W ~ —hjth of the mean free path of an upstream par¬ 
ticle. For this width the shock can still be considered 
almost step-like @- 

We model the downstream turbulence as a super¬ 
position of Alfen waves propagating parallel and anti¬ 
parallel to the plasma flow at constant speed. The 
turbulence spectrum has a power-law form with spec¬ 
tral index q for wavenumbers above some inverse cor¬ 
relation length k 0 . For k < kg the wave intensity 
per logarithmic bandwidth is assumed to be equal to 
the background field intensity, i.e., I{k) = B^k -1 for 
k < kq. In this work we use two values for q: 2, pro¬ 
ducing rigidity independent mean free paths, and 5/3, 
being consistent with the Kolmogorov phenomenology 
of turbulence. 

The downstream wave intensities can be calculated 
from known upstream parameters |e.g..liil.l?j. Hl|. and 
regardless of the cross helicity of the upstream wave 
field (only parallel or anti-parallel waves, or both), 
there are always both wave modes present in the 
downstream region; in this work we assume the cross 
helicity to vanish in the upstream. The amplification 
of the waves at the shock depends on the strength of 
the magnetic field as well as on the form of the tur¬ 
bulence spectrum, as shown in Fig. ^ where the am¬ 
plification of the different wave modes is shown as a 
function of (quasi-Newtonian) Alfvenic Mach number 

M = UA,i/ui- (2) 

«A,1 = /3i,H ( 3 ) 

is the upstream Alfven proper speed with /?a,i = 
va,i/c. The waves are seen to propagate pre¬ 


dominantly backward for relatively low-Mach-number 
shocks; this is the case for both the relativistic Jsj 
and the nonrelativistic 0 shocks. This enables the 
scattering center compression ratio rk to grow larger 
than the gas compression ratio r and, thus, to cause 
significantly harder particle spectra compared to the 
predictions of theories relying on fluctuations frozen-in 
to plasma flow §J fs]. As the Mach number increases, 
the downstream wave intensities approach equiparti- 
tion at the ultra-relativistic limit. 

Electrons scatter off the magnetic fluctuations res¬ 
onantly. The scattering frequency of electrons with 
Lorentz factor 7 is determined (see ref. 0 for details) 
by the intensity of waves at the resonant wavenumber 


k 


res — 


fie 

V 


^e,0 

C V 7 2 ~ 1 ' 


( 4 ) 


where fl e = is the relativistic electron gyrofre- 

quency and fl e ,o = 777 is its non-relativistic counter¬ 
part. 

Scatterings are elastic in the wave frame and the 
existence of waves propagating in both directions at 
a gi ven position, thus, leads to stochastic acceleration 
fl5l| . Since the spectrum of waves is harder below 
k = ko, scattering at energies 

Up n . , 

7 > 70 = -rf- > 1 ( 5 ) 


becomes less efficient. This suggests that the electron 
acceleration efficiency should decrease at 7 > 70 . In¬ 
stead of trying to fix the value of ko, we use a constant 
value 70 = 10 6 , which is in agreement with observa¬ 
tions of maximum Lorentz factor of electrons in some 
AGN jets [3. 

In addition to scattering, the particles are also as¬ 
sumed to lose energy via the synchrotron emission. 


3. RESULTS & DISCUSSION 

Simulations were run separately for low, intermedi¬ 
ate and high Alfvenic-Mach-number shocks (M = 3, 
M = 10, and M = 1000, respectively - see Fig. ^ for 
the corresponding wave intensity ratios), and for four 
cases of downstream turbulence; for turbulence spec¬ 
tral index q = 2 and <7 = 5/3 with downstream wave 
field calculated using wave transmission analysis de¬ 
scribed earlier, and with the downstream forward and 
backward waves being in equipartition. The proper 
speed of the shock is set to u\ = 10 c in all simula¬ 
tions. 

Some of the results (including all of the case of M = 
1000 for which the effects were, expectedly, barely vis¬ 
ible, and all those for which the downstream turbu¬ 
lence was calculated using the Alfven wave transmis¬ 
sion model) have been omitted from this paper due to 
limited space; all of them are, however, available at 
http://www.astro.utu.fi/red/qshock.html 
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Figure 2: Steady-state energy distribution of shock-injected particles at different distances from the shock. Upper 
panels show contours of log {E^X) as function of energy and location, while the lower panels present the energy 
distributions as slices at certain locations in the downstream. On the left hand panel is a case with a shock with 
Alfvenic Mach number M = 3 and for turbulence corresponding to the Kolmogorov turbulence, and on the right panel 
the same for M = 10 and turbulence corresponding to that created by the particles themselves. Equipartition is 
assumed between forward and backward waves. 131. 


3.1. Particles Injected at the Shock 

In the first case we studied the effect of stochas¬ 
tic acceleration on particles that have been already 
accelerated at the shock. We injected the particles 
with relatively small initial energy into the simulation 
right behind the shock and allowed them to continue 
accelerating via the stochastic process. This kind of 
injection simulates the case of some already-energized 
downstream particles returning into the shock, but re¬ 
moves need of processing the time consuming bulk of 
non-accelerating thermal particles. The high-energy 
part of the particle energy distribution - which we 
are interested in in this study - is similar, regardless 
of the injection energy. 

In the case of high Alfvenic Mach number (M = 
1000, corresponding to magnetic field Bq ~ 1.4 mG in 
a hydrogen plasma and ~ 46 p,G in a pair plasma) the 
contribution of the stochastic process to the energy 
distribution of the particles was, expectedly, very in¬ 
significant compared to that of the first-order acceler¬ 
ation at the shock. This was the case regardless of the 


applied turbulence properties, and because the wave 
intensities calculated from the Alfven wave transmis¬ 
sion analysis are very close to equipartition for high- 
Mach-number shock (see Fig. [Q, the difference to 
the case of explicitely assumed equipartition was, as 
expected, minimal. 

For shocks with Alfvenic Mach number sufficiently 
low (M = 10 and M = 3, corresponding to « 0.1 
0.5 G in a hydrogen plasma, and to ss 5-15 mG in 
a pair plasma) the effect was, on the contrary to the 
high-M case, very pronounced. The stochastic process 
began to re-accelerate particles immediately after the 
shock front, and the whole spectrum started to slowly 
shift to higher energies in all kinds of turbulence mod¬ 
els studied. Two examples for the case assuming the 
downstream wave intensities to be in equipartition are 
shown in Fig. [21 

When comparing otherwise similar cases that differ 
only for the downstream cross helicity (i.e., whether 
the wave field is resulting from the wave transmis¬ 
sion calculations [|j or an equipartition of parallel and 
anti-parallel waves is assumed), the calculated wave- 
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Figure 3: Same as Fig. Qbut for particles injected uniformly throughout the downstream region. 


El- 


transmission cases with more anti-parallel waves (see 
6 ] and FigQJ show stronger first-order acceleration, 
but weaker stochastic acceleration. This is because 
of the larger scattering center compression ratio in 
the wave-transmission case leading to more efficient 
first-order acceleration |J and Virtanen & Vainio, in 
preparation] and, on the other hand, faster momen¬ 
tum diffusion rate in the equipartition case leading to 
more efficient stochastic acceleration. 

3.2. Particles Injected Across the 
Downstream 

Next we assumed that a constant injection mecha¬ 
nism exists throughout the downstream region - this 
mimics a case where turbulent fluctuations cascade to 
higher wavenumbers and inject a fraction of the ther¬ 
mal electrons to the stochastic acceleration process. 
We injected particles at a constant initial energy uni¬ 
formly and isotropically within the whole downstream 
region. Examples of results are shown in Fig.j^for pa¬ 
rameters otherwise equal to those in Fig. |3 

The general behavior of the acceleration process is 
very similar for both injection methods: the piling up 
of the particles with 7 > 70 is visible and the effect of 


the turbulence spectrum is the same. For the constant 
injection, however, particles begin to form a power- 
law plateau in the energy range between the injection 
energy and 70 . The produced spectral index depends 
on the spectral index_of the magnetic field fluctuations 
q as a = q — 1 EMI; for (? = 2 , the particle spectral 
index er ~ 1, and for q = 5/3, a ~ 0.6. 

Also here the composition of the downstream wave 
field affects the resulting spectrum: in the case of 
downstream turbulence calculated using the wave 
transmission model the particle population immedi¬ 
ately behind the shock front extends to slightly higher 
energies than in the equipartition case, but for the 
latter the stochastic acceleration is clearly quicker. 
The former effect is due to the scattering center com¬ 
pression ratio being larger for lower -M shocks, while 
the latter is due to larger velocity differences in the 
equipartition wave field. 

3.3. ”Re-acceleration of Re-accelerated 
Particles” 

The prolonging of the mean free path of particles 
with 7 > 70 (or k res < ko) has, in addition to the 
drop of the stochastic acceleration efficiency, another 
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Figure 4: Stochastically accelerated particles return to the shock and get re-accelerated by the first-order process. On 
the right the particles are gathered at downstream free-escape boundary. 0 


interesting feature: particles that are already ener¬ 
gized first in the shock by the first-order mechanism 
and then in the downstream by the second-order ac¬ 
celeration, become able to return back to the shock 
and get ”re-injected” into the first-order acceleration 
process. Injection energies for this second injection 
are very high, or course, so particles get accelerated 
to ultra-high energies; this is seen as a distortion of 
the high-energy part of the left-hand contour panels 
in Figs. |3 and 0 and especially clear in Fig 01 where 
the low Mach number M = 3 and turbulence corre¬ 
sponding to q = 2 , both favorable for fast stochastic 
acceleration, allow particles to be accelerated to en¬ 
ergies of the order of 70 while still being sufficiently 
close to the shock front in order to being able to return 
to the first-order acceleration process. As a result, a 
clearly visible high-energy (quasi) power-law is seen 
in the upper part of the contour as well in the energy 
spectrum collected at the downstream border. 

3.4. Combination of Injection 
Mechanisms 

We also investigated what kind of particle energy 
spectra the two injection schemes - one operating at 
the shock, and another operating uniformly through¬ 
out the downstream region - are able to create. Fig.|3 
gives an example of combination of these. In simula¬ 
tions these two acceleration cases were mostly kept 
separate for the sake of simplicity, but there should 
be no reason to assume the separation be present also 
in nature. Also the relative amounts of shock- and 
downstream-injected particles are not fixed by our 
model, but instead considered a parameter. 


4. CONCLUSIONS 

We have studied stochastic particle acceleration in 
the downstream region of a relativistic parallel shock. 



Figure 5: Example of a combination (solid line) energy 
spectra of particles injected to the acceleration process at 
the shock (dotted) and throughout the downstream 
region (dashed). (Model parameters are q = 5/3 and 
M = 10.) The particles are collected at the downstream 
free escape boundary, ~ 4 x 10 12 cm away from the 
shock, and the number of the particles injected at the 
shock is larger than the number of uniformly injected 
particles by a factor of 100. 0 


Applying the wave transmission calculations of [§j and 
assuming the cross helicity to vanish in the upstream, 
we have modeled the turbulence of the downstream 
region as a superposition of Alfven waves propagat¬ 
ing parallel and anti-parallel to the plasma flow. Us¬ 
ing a kinetic Monte Carlo simulation we have modeled 
the second-order Fermi acceleration of electrons in the 
shock environment, and considered cases of accelera¬ 
tion of downstream-injected particles, as well as that 
of particles injected at the shock. We have shown that 
the stochastic acceleration can, indeed, have remark¬ 
able effects in both cases. This result is even more 
pronounced if the two downstream Alfven wave fields 
are assumed to be in equipartition. 

The behavior of the particle energy distribution in 
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the stochastic process depends heavily on the strength 
of the background magnetic field; in the cases of weak 
magnetic field and large quasi-Newtonian Alfvenic 
Mach number the effects of stochastic acceleration are 
faded out by the much stronger first order accelera¬ 
tion. Also the magnetic field turbulence spectrum af¬ 
fects the acceleration efficiency: for Kolmogorov tur¬ 
bulence with q = 5/3 the spatial scales are up to an 
order of magnitude shorter than in the case of q = 2 
turbulence. Although the spatial scales in simulations 
presented here are enormous compared to those asso¬ 
ciated with shock acceleration (the first-order process 
in the immediate vicinity of the shock front), in case 
of blazars and other active galactic nuclei the scales 
are still orders of magnitude too small to be resolved 
even in the VLBI pictures - regardless of the turbu¬ 
lence and used magnetic field strength. Also the ac¬ 
celeration time scales are short: the time required to 
shift the whole spectrum from the initial energy range 
to 7 buik > 10 6 ranged from 10 to 50 minutes in the 
M = 10 case, and for M = 3 the times were < 1 
minute, as measured in the shock frame. 

In the cases where the stochastic acceleration was 
quick enough for particles to reach the 70 energy while 
being still sufficiently close to the shock in order to 
be able to make their way back to the upstream re¬ 
gion due to their prolonged mean free path, the first- 
order mechanism was able to re-accelerate the return¬ 
ing high-energy particles to even higher energies. This 
led to forming of a new (quasi-)power-law at energies 
7 > 10 7 in some cases. 

One notable feature of the present model is that in 
the case of a uniform injection process in the down¬ 
stream region, power-law spectra with high and low 
energy cut-offs are formed. Depending on the turbu¬ 
lence, particle energy spectra have power-law spectral 
indices of 0.5-1 with lower and higher energy cut-offs 
at 71 S3 10 1 ~ injection and 72 ~ 10 6 = 70 , respec¬ 
tively. These particles would produce synchrotron 
spectra with photon spectral indices —0.5 < a < 0 
in the GHz-THz regime for various initial parame¬ 
ters. These properties are quite similar to those of 
flat-spectrum sources, for which typical spectra with 
a > —0.5 in the GHz region and flare spectra with 
a S 3 — 0.2 in the optically thin region of the spectrum 
are seen; see e.g. [lJj and references therein. 

To conclude, the main results of this study are: 

• Stochastic acceleration can be a very efficient 
mechanism in the downstream region of parallel 
relativistic shocks, provided that the magnetic 
field strength is large enough in order to make 
the Alfvenic Mach number approach the critical 
Mach number (M c = y/r) of the shock, i.e., to 
increase the downstream Alfven speed enough to 
allow for sufficient difference in speeds of parallel 
and anti-parallel Alfven waves required for rapid 
stochastic acceleration. 


• In the case of a continuous injection mechanism 
in the downstream region particle energy distri¬ 
butions with hard spectral indices can be formed 
between the injection energy and 70 . The accel¬ 
erated particle populations could produce syn¬ 
chrotron spectra very similar to those of flat 
spectrum sources. 

• The interplay between the first- and second- 
order Fermi acceleration at relativistic shocks 
can produce a variety of spectral forms not lim¬ 
ited to single power laws. 

Remark: We wand to emphasize that these simu¬ 
lations are based on the test-particle approximation, 
i.e., the effects of the particles on the turbulent wave 
spectrum and on the shock structure are neglected. 
Including these effects in a self-consistent manner may 
lead to notable effects on the resulting spectrum, see 
e-g. 0 . Including these effects to our model is, how¬ 
ever, beyond of the scope of the present simulations. 
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